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АЛГОРИТМ РАСЧЕТА РАВНОВЕСНЫХ СОСТАВОВ И МАСС ФАЗ 
В МНОГОКОМПОНЕНТНЫХ РАСТВОРАХ ЭЛЕКТРОЛИТОВ 

Постановка задачи. Расчет составов и масс раствора электролита и 
твердых фаз (постоянного состава) является важной технологической за¬ 
дачей в областях галургической промышленности, геохимии и гидрогео¬ 
логии. Эта задача тесно связана с расчетом диаграмм растворимости со¬ 
ответствующих систем, однако имеет специфику, поскольку диаграммы 
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растворимости определяют составы, но не массы сосуществующих фаз. 

Сформулируем общую задачу следующим образом. Пусть для системы 
в изотермо-изобарических условиях ( Т , Р=сопаІ) заданы числа молей 
компонентов (щ), а также набор всех возможных твердых фаз постоянно¬ 
го состава(1.. .А). Необходимо найти фазовый состав и массы фаз гетеро¬ 
генной системы, отвечающие состоянию устойчивого равновесия. 

Для решения такой задачи естественно воспользоваться критерием 
минимизации свободной энергии Гиббса гетерогенной системы — С» 
(см., например, работы [1-3]): 

(& 8 ) т ,Р, П) .-^тіп. (1) 

Суть применяющегося алгоритма такова. Пусть (в рамках той или иной 
модели) мы знаем зависимость химических потенциалов компонентов в 
растворе Нг. (г ) от переменных состава, а также химические потенциалы 
соединений в твердой фазе (щ (8) =щ (8)0 ). Тогда с использованием уравне¬ 
нии материального баланса всегда можно определить зависп- 

М °Г'І/ м ^ м \° Т масс твердых фаз гетерогенной системы: О е = 
=С 8 (Д/ 1 .. .ін*). Таким образом, задача (1) сводится к решению (к) урав¬ 
нений, нелинейных относительно М { . Если к достаточно велико (&>20), 
то решение становится весьма трудоемким и требует соответствующих 
характеристик ЭВМ (прежде всего быстродействия). Наличие граничных 
условий 0 также усложняет задачу. В настоящей работе излагается 
алгоритм расчета, использующий решение минимального числа уравнений. 

Термодинамические соотношения. Из термодинамики гетерогенных 
систем известно [4-6], что критерий равновесия гетерогенной системы 
(У, Р сопзі) эквивалентен системе уравнений фазового и химического 
равновесия (для фазы переменного состава раствора). Введем функции: 

Др,=щ ( ' ) ~щ (8)0 , (2) 

где щ <8)0 — химический потенциал і-го соединения в твердой фазе; щ (П - 
его химический потенциал в растворе, рассчитанный с помощью уравне- 
нии химического равновесия из химических потенциалов компонентов 
раствора. Тогда, если Др.сО, раствор ненасыщен; если Дщ>0 — пересы¬ 
щен; если Др,=0, раствор находится в равновесии с і-той твердой фазой. 

Таким образом, если гетерогенная система находится в состоянии 
устойчивого равновесия, то: 

для всех имеющихся твердых фаз (Л/*>0): Дщ=0: 

для всех отсутствующих твердых фаз (Л/,=0): Др-<0. 

Если хотя бы для одной фазы (к) Дщ>0, то система не находится в со¬ 
стоянии устойчивого равновесия и этой фазе выгодно образовываться из 
раствора согласно (1). 

Алгоритм расчета. Идея алгоритма основана на том, что потенциал 
Гиббса есть функция состояния и не зависит от того, каким путем мы 
придем к состоянию равновесия. Добавим в исходную систему содержа¬ 
щую щ молей компонентов, п<*> молей воды, так чтобы равновесная систе¬ 
ма отвечала гомогенному раствору. Другими словами, для всех к твердых 
фаз справедливо: 

ДМі<0. (3) 

1) Теперь малыми шагами Д п ш <п® будем выводить из системы п$ 
молей воды. На каждом этапе проверяется условие: 

Дц<> 0. (4) 

2) Пусть на некотором этане выведения для фазы а условие (4) ока¬ 
залось выполненным. Тогда последующее выведение будет происходить 
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яр* хемостатировании раствора по соединению а: 

ДЦа=0. (5) 

Таким образом, мы будем одновременно удалять из раствора воду и 
фазу ’ ®' ^ вигаться яо гиперобъему кристаллизации фазы а на диаг¬ 

рамме растворимости ( М а увеличивается). При этом условие (4) будет 
проверяться для всех остальных к— 1 твердых фаз. 

3) Пусть на некотором этапе условие (4) оказалось выполненным для 
фазы о. Іогда последующее выведение будет происходить при хемостати¬ 
ровании раствора по двум соединениям — а и Ъ: 

Др о =ДіХб=0. (6) 

Таким образом, мы будем удалять из раствора воду и фазу Ь ( М ь уве¬ 
личивается ). В это время фаза а будет либо также удаляться из раствора 
(М а увеличивается, эвтоническая кристаллизация), либо будет происхо¬ 
дить ее растворение ( М а уменьшается, перитоническая кристаллизация), 
т. е. мы будем двигаться по гиперобъему совместностной кристаллизации 
фаз а и Ъ. Если до выполнения условия (4) для третьей фазы с растворит¬ 
ся вся ранее образованная фаза а, то мы оставляем только одно хемоста- 
тирование: 

Дць=0, (?) 

т. е. переходим в более вариантный гиперобъем кристаллизации одной 
фазы Ь. 

4) Пусть этого не произошло, и на определенном этапе выведения усло¬ 
вия (4) оказалось выполненным для фазы с. Тогда дальнейшее хемоста- 
тирование пойдет сразу по трем фазам: 

Др, а ==Дрь=Д|х с =0, (8) 

что соответствует решению системы трех нелинейных уравнений фазового 
и химического равновесия. Вода и фаза с будут Удаляться из раствора 
(М с увеличивается), фазы а и Ъ будут либо совместно удаляться из рас¬ 
твора (М а Мь увеличиваются, эвтоническая кристаллизация), либо одна 
из них будет удаляться, а другая растворяться (М а и изменяются анти- 
батно перитоническая кристаллизация I рода), либо обе они будут рас¬ 
творяйся Тма, М ь уменьшаются, перитоническая кристаллизация II рода). 
Полному растворению одной на фаз будет отвечать переход . более ва- 

риантный гиперобъем: 

Дц с = Д|Х ь ==0 или Д|х <: =Ара=0, (9) 

_ ѵ <Ьяя — пеоеход в еще более вариантный ги- 
а полному растворению обеих фаз переход в ощ 

перобъем: / 4ПЧ 

Др с =0. 

5) Процесс удаления из раствора продолжается таким же способом до 

полного выведения моЛе * следующими свойствами: 

Приведенный алгоритм <^ ла Д ся устойчивое состояние равнове- 

в Результате поиска всегд Д ой ерки решения на устои- 

си* гетерогенной системы и дополни ^ делать п0 алгоритмам, 

чивость проводить не нужно (как э 

изложенным в работах [1,3]); „„ лИ „„ ргкие Т ак и перитонические рав- 
позволяет рассчитывать как э сразу с двумя различными 

яовесия любого рода, равновесия Р аств °Р° ва И Р Т „, В алгоритме учтено и 
кристаллогидратами одного оолево фазам, и «расхемостатирова- 

хемостатирование по образованным твердым ц, ш 


ние» в случае растворения последних, причем процесс хемостатирования - 
«расхемостатирования» нерекурсивен в отличие от алгоритма, приведен¬ 
ного в работе [1]! 

не допускает (в отличие от работы [ 1 ]) неопределенности с порядком 
хемостатирования и количеством хемостатируемых фаз. В приведенном 
алгоритме новые твердые фазы появляются или исчезают строго пооче¬ 
редно в соответствии с движением по фазовой диаграмме при изотермиче¬ 
ском удалении воды из раствора. 

Таким образом, никогда (даже формально) не может быть нарушено 
правило фаз Гиббса. 

Если бы в качестве начального состояния гетерогенной системы мы 
выбрали неразбавленный гомогенный пересыщенный но неопределенному 
числу твердых фаз раствор (см., например [ 1] Т то алгоритм не обладал 
бы приведенными свойствами. Действительно, пусть такой раствор ока¬ 
зался пересыщенным по р твердым фазам, а максимальное число равно¬ 
весно сосуществующих с раствором твердых фаз равно д<р- Тогда в за¬ 
висимости от выбора алгоритм придет к одному из С р ч решений со всеми 
признаками устойчивого равновесия, в то время как реальное устойчивое 
равновесие может не отвечать ни одному из этих решений. 

Настоящий алгоритм реализован нами для пятикомпонентной (взаим 
ной) галургической системы !Ча + , М§ 2+ , Са 2+ ||С1“, 80 4 2_ -Н 2 0 для интер¬ 
вала температур — 20-з+30° С. Коэффициенты активности задавались е 
использованием уравнений Питцера [7]. обеспечение которых осущест¬ 
влялось циклом подпрограмм типа [8]. Точность расчетов по приведен¬ 
ному алгоритму, естественно, не отличается от точности расчета диаграмм 
фазовых равновесий, широкий класс которых весьма точно прогнозируется 
по методу Питцера [9-13]. 
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